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The  realistic  evsduation  of  weapon  system  performance  is  a  critical  element 
in  DoD  decision  making.  Since  many  sensors  and  systems  cannot  operate  in 
clouds,  the  assessment  of  cloud  limitations  on  their  performance  plays  a  key  role 
in  the  evaluation  process.  Therefore,  the  DoD  needs  cloud  models  to  support  the 
operational  evaluation,  acquisition,  studies  and  analysis,  and  war  gaming  com¬ 
munities.  The  development  of  reliable  cloud  models  requires  accurate  cloud  cli¬ 
matologies.  The  purpose  of  this  atlas  is  to  establish  a  global  cloud  climatology 
that  details  seasonal  and  diurnal  variations. 

Since  cloud  climatologies  are  critical  elements  in  global  climate  modeling, 
monitoring  of  climatic  variability,  and  earth-oriented  space  experiments,  a  broad 
spectrum  of  the  non-DoD  meteorological  community  has  been  interested  in  this 
topic  in  recent  years.  The  establishment  of  the  International  SateUite  Cloud  Cli¬ 
matology  Project  (ISCCP)  indicates  wide  recognition  of  the  need  for  an  accurate 
global  cloud  climatology.  According  to  Schiffer  and  Rossow,  ^  the  goal  of  ISCCP 
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is  to  estimate  the  global  distribution  of  cloudiness  by  collecting  and  analyzing 
satellite  data. 


2.  BACKGROUND 
2.1  Previous  Work 

A  variety  of  global  and  regional  atlases  have  already  been  compiled  using 

2 

numerous  data  sources.  The  Handbook  of  Meteorology  contains  global  analyses 

0 

of  monthly  mean  cloudiness  for  every  other  month  starting  with  January.  London 
compiled  a  mean  seasonal  global  cloud  climatology  for  the  Northern  Hemisphere 

4 

based  on  surface  observations  taken  at  1200  LST.  Sadler  used  satellite  data  to 

5 

estimate  the  average  monthly  cloud -cover  for  the  tropics.  Sherr  et  al  compiled 
a  master  file  of  tabulated  cloud  statistics  and  cloud  distributions  on  a  world-wide 
basis,  including  diurnal  and  seasonal  changes.  They  defined  29  large-scale  ho¬ 
mogeneous  cloud  regions  such  that  the  statistics  from  a  single  station  within  the 
region  would  adequately  represent  the  entire  region.  Satellite  imagery  was  in- 

C 

eluded  among  the  data  sources.  Hall  analyzed  cloud  characteristics  at  32  loca¬ 
tions  in  Russia.  Frequencies  of  coverage  were  given  by  cloud  types  -  low,  mid- 

7 

die,  and  high.  Miller  and  Feddes  compiled  monthly,  seasonal,  semiannual,  and 
annual  global  charts  of  relative  cloud-cover  using  meteorological  satellite  data 

g 

acquired  during  the  period  1  January  1967  through  31  December  1970.  Malberg 


2.  Berry,  F.  A.,  Jr.,  Bollay,  E.,  and  Beers,  N.  R.  (1945)  Handbook  of  Mete¬ 

orology,  McGraw-Hill  Book  Company,  Inc.,  New  York,  pp.  988-9^6. 

3.  London,  J.  (1957)  A  Study  of  the  Atmospheric  Heat  Balance,  AFCRC-TR-57- 

287,  Final  Report,  Contract  No.  AF  l6(l22)-165.  New  York  University, 
ASTIN  No.  11  7227. 

4.  Sadler,  J.  C,  (1966)  Average  Monthly  Cloud  Cover  for  the  Global  Tropics  as 

Determined  from  Satellite  Observations.  Preprint  Pamphlet,  Contract  iJo. 
AF  19  ^626)-3860,  Hawaii  Institute  of  Geophysics,  University  of  Hawaii. 

5.  Sherr,  P.  E.,  Glasen,  A.  H. ,  Barnes,  J.  C. ,  and  Willand,  J.  H.  (1968) 

World-Wide  Cloud  Cover  Distribution  for  Use  in  Computer  Simulations, 
NASA  CR-61226,  Final  Report,  Contract  No.  NAS  8-21040,  Allied  Re  - 
search  Associates,  Inc. ,  Concord,  Mass. 

6.  Hall,  R.  J.  (1969)  Cloud  Cover  of  the  USSR.  SAMSO-TR-69-147,  Vol.  1, 

Contract  No.  F04701-68-C-0200,  Aerospace  Corp. 

7.  Miller,  D.  B. ,  and  Feddes,  R.  G.  (1971)  Global  Atlas  of  Relative  Cloud 

Cover  1967-70.  NOAA/NESS.  USAFETAC. 

8.  Malberg,  H.  (1973)  Comparison  of  mean  cloud  cover  obtained  by  satellite 

photographs  and  ground-based  observations  over  Europe  and  the  Atlantic, 
Mon.  Wea.  Rev.  101:893-897. 


compared  satellite  and  ground-based  observations  of  mean  cloud-cover  for  Europe 

g 

and  the  Atlantic.  Berlyand  and  Strokina  present  a  world  atlas  of  mean  monthly 
cloud  amounts  for  every  other  month  starting  with  January,  with  data  taken  from 
approximately  3400  surface  stations  and  from  satellite  observations.  The  maps 
represent  cloud-cover  distributions  for  daylight  hours  only.  Hahn  et  al^*^  pub¬ 
lished  an  atlas  of  co-occurrence  of  different  cloud  types  over  the  ocean.  Ship  ob¬ 
servations  were  used,  and  the  clouds  grouped  into  six  types.  The  results  are  ex¬ 
pressed  as  contingency  probabilities  at  a  resolution  of  15“  latitude  X  30“  longitude. 
They  also  published  a  subsequent  atlas on  the  co-occurrence  of  cloud  types  over 

land  using  a  similar  analysis  procedure  on  land  observations.  Henderson-Sellers 
12 

and  Hughes  produced  a  one-year  (1979)  global  cloud  climatology  using  3DNEPH 
data  exclusively. 

Several  authors  have  modeled  the  frequency  distribution  of  cloudiness  as  seen 
by  the  ground  observer  or  satellite.  Each  model  includes  all  fractions  of  cover¬ 
age  ranging  from  clear  to  overcast,  and  is  claimed  to  have  versatile  statistical 

characteristics  to  simulate  J-  or  U-shaped  distributions  of  sky  cover  as  well  as 

13 

more  conventional  bell -shaped  distributions.  Henderson-Sellers  et  al,  in  ex¬ 
plaining  the  need  for  a  global-scale  cloud  archive,  recommend  the  Beta  distribu- 

14  15 

tion  model  as  developed  by  Falls,  '  Its  density  function  is 

f(x)  =  x^'^  (1  -  a,b>  0  ,  (1) 


9.  Berlyand,  T.  G.  ,  and  Strokina,  L.  A.  (1980)  Global  distribution  of  total 
cloud  amount,  Gidrometeoizdata,  Leningrad  (translated  by  S.  Warren). 

10.  Hahn,  C.  J. ,  Warren,  S.  G. ,  London,  J. ,  Chervin,  R.  M. ,  and  Jenne,  R. 
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Prediction  Division,  Boulder,  Colo. 

11.  Hahn,  C.  J. ,  Warren,  S.  G. ,  London,  J. ,  Chervin,  R.  M.  ,  and  Jenne,  R. 

(1984)  Atlas  of  Simultaneous  Occurrence  of  Different  Cloud  Types  Over 
Land,  NCAR  Technical  Note  No.  TN-241+STR,  Atmospheric  Analysis  and 
Prediction  Division,  Boulder,  Colo. 

12.  Henderson-Sellers,  A.,  and  Hughes.  N.  A.  (1984)  Cloud  Archiving  Strate¬ 

gies,  AFGL-TR-84-0166,  Final  Report,  Contract/Grant  No.  AFOSR-83- 
0118,  University  of  Liverpool,  Liverpool,  U.  K. 

13.  Henderson-Sellers,  A.,  Hughes,  N.  A.,  and  Wilson,  M.  (1981)  Cloud  archiv¬ 

ing  on  a  global  scale;  a  discussion  of  principles.  Bull.  Am.  Meteorol.  Soc. 
62:1300-1307. 

14.  Falls,  L.  W.  (1973)  The  Beta  Distribution:  A  Statistical  Model  for  World 

Cloud  Cover.  NASA  Technical  Memorandum  No.  64714,  Marshall  Space 
Flight  Center,  Ala, 

15.  Falls,  L.  W.  (1974)  The  beta  distribution:  a  statistical  model  for  world  cloud 
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where  r  is  the  gamma  function,  x  is  the  fraction  of  sky  covered,  and  a  and  b  are 
parameters  whose  values  can  vary  to  permit  varying  characteristics  of  the  cloud - 
cover  distribution.  Falls  assigned  values  to  a  and  b  in  29  regions  of  the  world. 
While  the  density  function  is  easy  to  compute,  the  cumulative  probability  is  diffi¬ 
cult  to  obtain  since  no  succinct  analytical  expression  exists.  Therefore,  tables, 
numerical  integration,  or  other  approximate  methods  are  necessary  to  obtain 
cumulative  probabilities. 

1 6 

Bean  and  Somerville  also  used  the  Beta  distribution  to  model  cloud-cover 
distributions  world-wide.  They  used  the  infrared  measurements  from  a  TIROS 
satellite  taken  over  a  period  of  approximately  four  years  and  determined  the  pro¬ 
portion  of  cloud  cover  world-wide  over  2,  5*  latitude -longitude  grids.  Estimates 
of  a  and  b  were  then  obtained  and  grouped  into  homogeneous  regions.  Thus,  the 

globe  was  divided  into  homogeneous  cloud -cover  regions  for  each  season. 

17 

Somerville  et  al  used  the  Johnson  Sg  family  of  distributions  to  model  total 
sky-cover  for  seven  diverse  locations.  The  Sg  family  is  given  by 

Z  =  y  +  T7  in  [x/(l  -  x)]  .  (2) 

The  variable  x  is  equal  to  the  fractional  coverage,  and  the  parameters  y  and  p  are 
determined  by  the  data.  The  variable  Z  is  the  equivalent  normal  deviate  (END)  of 
the  cumulative  probability  of  x.  Z  has  a  mean  of  zero  and  a  standard  deviation  of 
1.  The  solution  gives  the  probability  that  the  sky  cover  will  be  less  than  or  equal 
to  the  fractional  cover  x.  Someiwille  et  al  developed  96  pairs  of  values  for  y  and 
n  for  each  station  (8  times  per  day,  12  months)  using  simple  linear  regression. 

The  root  mean  square  (rms)  error  for  all  stations  ranged  from  0.  0175  to  0.  0595. 

Somerville  and  Bean  developed  another  cloud  model  using  the  S-distribution, 
which  is  given  by 

F(x)  =  1  -  (1  -  x®)^  o.jS  >  0  .  (3) 

As  in  Eqs.  (1)  and  (2),  x  is  equal  to  the  fractional  coverage  and  F(x)  is  the  proba¬ 
bility  that  the  sky  cover  is  less  than  x.  The  parameters  a  and  0  can  be  estimated 
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iteratively  using  either  the  method  of  maximum  likelihood  or  a  nonlinear  regres¬ 
sion  technique.  Both  methods  are  considerably  more  complicated  than  the  simple 
linear  regression  used  to  obtain  Eq.  (2).  Somerville  developed  96  different  pairs 
of  values  for  a  and  3  for  23  stations  around  the  world.  The  rms  errors  ranged 
from  0.003  to  0.071.  Strictly  speaking,  this  model,  as  well  as  the  previously 
described  models,  yields  zero  for  the  probability  of  all  clear  and  full  overcast; 
hence  they  cannot  be  theoretically  accurate.  The  difficulty  is  overcome,  in  prac¬ 
tice,  by  stipulating  that  the  integration  for  clear  be  taken  from  0  to  0.  5  tenths,  and 
for  overcast  from  9.  5  to  10  tenths. 

2.2  Current  Work 

This  AFGL  atlas  uses  a  mathematical  model  that  describes  the  frequency  dis¬ 
tribution  of  total  sky-cover.  The  model,  known  as  the  Boehm  Sawtooth  Wave 
(BSW)  simulation  model,  is  explained  in  more  detail  in  Section  3.  Like  the  models 
of  Falls  and  Somerville,  the  BSW  model  uses  two  parameters.  Our  model  has 
several  advantages  over  theirs,  however.  First,  both  parameters  have  physical 
meaning.  One  parameter  is  the  mean  sky-cover,  which  we  identify  as  the  single - 
point  probability  of  a  cloud.  The  other  parameter,  which  we  call  scale  distance, 
is  the  distance  over  which  there  is  a  correlation  coefficient  of  0.  99.  The  corre¬ 
lation  coefficient  between  points  separated  by  greater  distances  can  be  found  by 
Eq.  (4)  (Section  3).  It  has  been  found  that  the  scale  distance  varies  systematically 
with  time  of  day  and  time  of  year,  which  enhances  its  choice  as  a  parameter.  The 
second  advantage  of  the  BSW  model  is  that  it  can  be  extended  to  provide  the  distri¬ 
bution  of  cloud  cover,  from  clear  to  overcast,  over  areas  of  varying  sizes,  from  a 
2  2 

few  km  to  100, 000  km  .  A  third  advantage  is  that  conditions  of  absolute  clear 
and  full  overcast  have  meaningful  probabilities,  estimated  unequivocally  by  the 
model. 

The  scale  distance,  in  combination  with  the  mean  sky-cover,  permits  the  cal¬ 
culation  of  the  sky-cover  distribution  from  clear  to  overcast.  Thus,  the  frequency 
distribution  of  sky  cover  can  be  found  for  locations  world-wide  using  the  global 
analyses  of  mean  sky-cover  and  scale  distance  given  in  the  atlas.  Twenty  maps 
are  presented  for  each  of  these  parameters,  one  set  four  times  daily  (0000-0200, 
0600-0800,  1200-1400,  and  1800-2000  LST)  for  each  mid-season  month  (January, 
April,  July,  and  October),  and  another  set  of  monthly  means  for  each  mid-season 
month.  The  model  provides  for  convenient  data  storage,  since  two  numbers  define 
the  entire  frequency  distribution. 
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3.  THE  BSW  MODEL 


The  Boehm  Sawtooth  Wave  (BSW)  simulation  model  was  developed  by  Maj. 

Albert  R.  Boehm  of  AFGL.  The  model,  described  in  detail  in  Burger  and 
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Gringorten,  generates  fields  (or  maps)  of  normally  distributed  values  called 
equivalent  normal  deviates  (ENDs).  An  analytical  solution  has  been  found  for 
BSW  for  the  spatial  correlation  function  p(s)  between  points  separated  by  distance 
(s),  in  units  of  the  wavelength; 

p(s)  =  1  -  (12/if)s  +  3s^  for  0  ^  s  s  1  (4) 

=  1  -  (12/ir)s  +  3s^  +  (24/ir)[cos'^(l/s)  -  -  1)  forl<s:<2. 

where  s  is  the  standardized  distance  given  by 

s  =  s'/A  ,  (5) 

s'  is  the  measured  distance  (km),  and  A  is  the  wavelength  (km).  For  s  >  2,  p(s) 
oscillates  between  negative  and  positive  values,  converging  on  zero.  The  graph 
of  p(s)  vs  s  is  shown  in  Figure  1.  The  correlation  is  high  for  short  distances,  but 
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DISTANCE  (UNITS  OF  WAVELENGTH) 


Figure  1.  Distance  vs  Correlation  Coefficient 
for  BSW  Model 


drops  monotonically  to  zero  at  s  =  0.3684,  becomes  negative,  and  crosses  zero 
again  at  s  °  0.  9049.  The  distance  over  which  p(s)  decreases  to  0.  99  is  the  model 
parameter,  scale  distance  (r/km).  Eq.  (4)  yields  p(s)  ^  0.  99  over  a  distance  of 


19.  Burger,  C.  F.,  and  Gringorten,  I.  I.  (1984)  Two-Dimensional  M^eling  for 
Lineal  and  Areal  Probabilities  of  Weather  Conditions.  AFGL-TR-84-0126. 
AD  A147970. 


1  km,  (r  =  1  km)  if  A=  340  km.  The  ratio  between  r  and  A  is  constant,  chosen 
for  this  atlas  to  be  1/340.  Thus, 

r  =  (A/ 340)  km  .  (6) 

19 

As  previously  described,  an  analytical  solution  for  the  probability  of  a  max¬ 
imum  condition  in  an  area  or  fraction  of  an  area  has  not  yet  been  found.  Instead, 
we  have  obtained  approximate  solutions  by  a  Monte  Carlo  method.  The  results 
were  put  in  graphical  form,  one  graph  for  each  tenth  of  coverage,  from  0  to  10 
tenths.  The  graphs  for  0  through  5  tenths  are  shown  in  Figures  2  through  7.  The 
graphs  for  6  through  10  tenths  are  not  shown  because  they  are  merely  the  inverted 
mirror  images  of  the  graphs  for  4  through  0  tenths,  respectively.  On  the  graphs, 
the  vertical  axes  are  uniform  in  the  END  (y)  of  the  cumulative  probability  (P). 

The  horizontal  axis  has  a  uniform  z -scale,  where 


from  which  we  obtain 


A  is  the  area  of  interest. 


The  graphs  were  fitted  by  empirical  equations,  making  it  possible  to  obtain 
solutions  quickly  by  computer.  The  solution  is  as  follows; 

y'  =  a  +  ^  +  7  +  6  ,  (9) 

where  y'  is  the  END  of  the  cumulative  probability  of  sky  cover  for  a  given  frac¬ 
tional  coverage  (F/lO);  y'  will  vary  with  r  and  mean  sky-cover  P^.  The  terms  a, 
P,  7,  and  6  are  given  in  Table  1,  and  the  associated  constants  in  Tables  2  and  3. 
The  variable  y^  (from  Table  1)  is  the  END  of  1  -  P^j.  When  F  >  5,  the  signs  of  y' 
and  y^  are  reversed,  and  F  =  10  -  F.  In  addition,  when  z>  7  and  F  is  not  an  in¬ 
teger,  y'  -  y'  [INT(F)1  +  [F  -  INT(F)1  {y'[INT(F)  +  1]  -  y'[INT(F))  }.*  These 
solutions  are  valid  only  within  the  bounds  of  the  graphs,  and  therefore  include  the 
following  restrictions:  |y'  |  ^  3.72,  |y^|  <  4.7,  and  -1  <  z  <  8.  Note  that  the 
formula  permits  solutions  for  any  value  of  F  instead  of  just  the  integer  values,  as 
with  the  graphs.  In  other  words,  the  empirical  solution  interpolates  between 
graphs,  so  that  it  is  possible,  for  example,  to  use  the  model  for  eighths  of  cloud 
cover  as  well  as  tenths. 


"iNT  is  a  function  that  acts  to  truncate  any  number  to  its  integer  value. 


j.'  •  j.'*" 


7 


Figure  2.  Graph  for  Estimating  P{0/ 10)  and  P(<  10/ 10)  Coverage.  Each 
curve  corresponds  to  a  single  yg  value,  and  z  =  In  (V  A/r)/ln  2.  The  addi¬ 
tional  solid  lines  and  dashed  curve  illustrate  an  example  from  Section  3.  1 
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Figure  3.  Graph  for  Estimating  P{s  l/lO)  and  P{s  Coverage.  Each 

curve  corresponds  to  a  single  y^,  value,  and  z  >=  In  (^/T/r)/ln  2 
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Figure  5.  Graph  for  Estimating  P{s  3/10)  and  P{s  '^0)  Coverage.  Each 
curve  corresponds  to  a  single  value,  and  z  =  ln(^/X/^)/ln  2 


r/r)/ln  2 


Figure  6.  Graph  for  Estimating  P(<4/10)  and  P(<6/10)  Coverage.  Each 
curve  corresponds  to  a  single  value,  and  z  =  In {^/A/^)/ln2 


Table  2.  Constants  for  Eq.  (9)  (With  Table  3) 


n 

'^n 

0 

0.0000 

0.0000 

0.00 

1 

-0.0071 

0. 0086 

0.  53 

2 

-0.0205 

0.0097 

0.  52 

3 

0.0229 

0.0045 

0.64 

4 

-0.0517 

0. 0224 

0.45 

5 

0.0260 

0.  0087 

0.  71 

Table  3.  Constants  for  Eq.  (9)  (With  Table  2) 


m  =  1 

m  =  2 

B 

m 

m 

g(F) 

®  m 

m 

m 

g(F) 

^  m 

0.  1420 

-2.4090 

-0.44 

0. 1908 

-3.  8406 

1.0686 

-2.2854 

-0.46 

1.  8710 

-4.2412 

1. 9440 

-2.7842 

-0.41 

3. 0540 

-4.4752 

-0.  39 

2.2560 

-2.7356 

-0.40 

3.3960 

-4. 5228 

-0.  50 

2.0328 

-2.2910 

-0.41 

-2.0690 

-2. 5702 

-0.75 

5 

1.2356 

-1.2300 

_ 

-1.30 

1.6124 

_ 

-1. 5434 

-1.  51 

4.  PROCEDURE 

We  now  apply  the  results  of  Section  3  to  the  atlas,  giving  a  step-by-step  pro¬ 
cedure  for  using  both  the  graphs  and  the  equations.  Initially,  we  make  two  impor¬ 
tant  assumptions.  First,  we  adopt  the  mean  sky-cover  as  the  probability  of  a 
cloud  presence  vertically  above  a  single  point  (Pp)-  The  no-cloud  probability  then 
becomes  P^'  =  1  -  P^.  Second,  we  assume  that  the  celestial  dome,  as  seen  by  the 

weather  observer,  has  a  radius  of  27.8  km  (15  nmi),  giving  a  floor  space  of  area 
2 

A  =  2424  km  .  With  these  two  assumptions,  we  can  calculate  cumulative  proba¬ 
bilities  of  fractional  sky-cover  as  seen  by  a  ground  observer  (in  other  words,  the 
climatology  of  total  sky-cover)  from  the  scale  distance  and  mean  sky-cover  values 
in  the  atlas. 


4.1  Graplu 

The  graphical  procedure  for  obtaining  the  cumulative  probabilities  of  frac¬ 
tional  sky-cover  is  shown  in  the  following  steps. 

1.  Obtain  r  and  P„  at  the  desired  location  and  time  from  the  atlas  maps. 

o  2 

2.  Calculate  P^'  =  1  -  P^,  and  calculate  z  using  Eq.  (7)  with  A  =  2424  km  . 

3.  For  F  =  0-  to  5 -tenths  coverage: 

a.  Choose  the  graph  corresponding  to  the  appropriate  coverage,  O/lO 
through  5/ 10  (Figures  2  through  7). 

b.  Find  PJ  on  the  left-hand  probability  scale.  This  defines  the  y^  (END) 
curve. 

c.  Follow  the  y  curve  until  it  intersects  the  z  value. 

o 

d.  Move  left  horizontally  to  the  vertical  axis  and  read  P(<  F/lO). 

4.  For  F  >  5 -tenths  coverage: 

a.  Choose  the  graph  corresponding  to  the  complement  of  the  appropriate 
coverage. 

b.  Find  the  complement  of  P^'  on  the  left-hand  probability  scale.  This 

defines  the  y  curve, 
o 

c.  Follow  the  y  curve  until  it  intersects  the  z  value. 

•'o 

d.  Move  right  horizontally  to  the  vertical  axis  and  read  P(sF/lO).  (In  the 
case  of  10  tenths,  we  have  P(<  10/ 10)  and  P(overcast)  =  1  -P(<  10/10). ) 

We  illustrate  the  graphical  method  with  a  specific  example.  Suppose  we  are 
interested  in  the  sky-cover  climatology  at  a  certain  location  for  January  noontime. 

Step  1.  We  obtain  r  and  P^  from  the  appropriate  spot  on  the  January  1200- 
1400  LST  maps  (Figures  15(a)  and  15(b)].  Let  us  say  r  =  2.23  and  P^  =  0.63. 

Step  2.  Next  we  calculate  the  no-cloud  probability  P^'  and  1  -  P^  =  0.  37,  and 
z  =4.46  using  Eq.  (7),  with  A  =  2424  km^. 

Step  3.  To  estimate  the  climatic  probability  of  clear  skies,  use  the  graph 
for  O/lO  coverage  (Figure  2).  On  the  left-hand  axis,  find  P^'  =  0.  37,  Follow  the 
appropriate  y^  curve  (in  this  case  y^  =  0.  33,  the  END  of  P^',  as  shown  by  a 
dashed  line  in  Figure  2)  to  its  intersection  with  z  =  4.46  (solid  vertical  line). 

From  that  poi.  move  left  horizontally  (solid  horizontal  line)  to  the  vertical  axis 
and  read  the  pr,  ^ability  of  clear  skies,  P(0/ 10)  =  0.  165. 

To  estimate  the  probability  of  s  I -tenth  sky-cover,  go  to  the  graph  for  l/lO 
coverage  (Figure  3).  Using  ihe  same  procedure.  Step  3,  with  P^'  =  0.  37  and 
z  =  4.46,  we  obtain  P(:£  1/ 10)  =  0.235.  Similarly,  the  probability  of  «  2 -,  s3-, 
<4-,  and  ^  5-tenths  sky-cover  is  estimated  from  the  graphs  for  2/10  through 
5/ 10  coverage  (Figures  4  through  7),  The  values  are  shown  in  Table  4,  row  1. 

Step  4.  Recall  that  the  graphs  for  0  through  4  tenths  are  the  inverted  mirror 
images  of  the  graphs  for  10  through  6  tenths,  respectively.  Therefore,  we  esti¬ 
mate  P(:£6/l0)  coverage  with  the  graph  for  1  -  6/10  =  4/10  coverage  (Figure  6). 
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Table  4.  Cumulative  Probabilities  for  Each  Tenth  of  Cloud  Cover  From  the 
Example  r  =  2.23  and  Pq  =  0.  37  by  (1)  the  Graphical  Method  and  (2)  by  Com¬ 
puter 

0  <1  <2  <3  <4  <5  <6  <7  <8  <9  <10 

1  0.165  0.235  0.270  0.310  0.330  0.370  0.410  0.430  0.460  0.510  0.620 

2  0.169  0.234  0.270  0.305  0.333  0.364  0.396  0.426  0.466  0.512  0.603 


The  no-cloud  probability  becomes  P^'  =  1  -  P^'  =  1  -  0.  37  =  0.  63,  and  therefore 
y^  =  -y^  =  0.  33.  Now  follow  the  curve  for  y^=  0.  33  in  Figure  6  to  its  intersec¬ 
tion  with  z  =  4.  46.  From  that  point,  move  right  horizontally  to  the  vertical  axis 
and  read  the  probability  of  6-tenths  coverage,  P(<6/10)  =  0.410.  This  is  the 
complement  of  the  value  that  would  be  read  by  moving  to  the  left-hand  scale. 

To  estimate  the  probability  of  7 -tenths  sky-cover,  use  the  graph  for 
1  -  7/10  =  3/10  coverage  (Figure  5).  Using  the  same  procedure.  Step  4,  with 
Pq'  =  0.  63  and  z  =  4.46,  we  obtain  P(:£7/l0)  =  0.430.  Similarly,  the  probability 
of  ^8-,  S9-,  and  <  10-tenths  sky-cover  is  estimated  from  the  graphs  for  2/10 
through  O/lO  coverage  (Figures  4  through  2),  and  are  shown  in  Table  4.  The  final 
value,  P(<  10/ 10)  coverage,  is  the  probability  of  all  fractions  of  coverage  except 
overcast.  In  this  case,  P(<  10/10)  =  0.620,  and  therefore  the  probability  of  over¬ 
cast  is  P(overcast)  =  1  -  0.620  =  0.380. 


4.2  By  Computer 

Continuing  with  this  example,  we  now  solve  using  the  empirical  equations. 
The  method  is  simple  to  program  in  BASIC  or  FORTRAN.  We  begin  with  the 
same  values  of  P^,  r,  and  sky-dome  area  A,  and  compute  P^'  =  0.37  and  z  =4.46 
as  before.  We  then  proceed  as  follows; 

1.  Calculate  the  END  (y^)  of  P^'  using  the  National  Bureau  of  Standards 
(NBS)  formula 


y.  =  k  t- 


^o  ^1^ 

1  +  bjt  +  bgt^ 


+  e 


where 

k  =  -  1,  t  =  4/ln  — when  P  '  S  0.  5  , 

1  {p^r  ° 

k  =  1.  t  =4/ In - 1 - 5  whenP  '>  0.5  , 

1  (1  -P  T  ° 

o 


(10) 
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a  =  2.30753  b,  =  0.99229 
o  1 

aj  =  0.27061  bg  =  0.04481  . 

2.  Next,  choose  the  fractional  coverage  <F/10)  of  interest.  F  can  be  any 
number  in  the  range  0  :£  F  rs  10. 

3.  Enter  y^,  F,  and  z  into  Eq.  (9).  The  expressions,  conditions,  and  con¬ 
stants  for  the  equation  are  shown  in  Tables  1  through  3,  as  explained  in  Section  3. 
The  result  is  the  END  of  P(  <  F/  10),  y'. 

4.  Compute  P(:sF/10)  coverage  from  y*  using  the  NBS  formula 

P(<F/10)  =  k+0.5m  (l  +  Cjjy'l  +  C2y’^  +  Cg  1  y'  1  ^  +  c^y''^  j  +  e  , 
where 

k  =  0,  m  =  1  for  y'  <  0  , 
k  =  1,  m  =  -1  for  y'  >  0  , 
le|  <  0.00025  , 

Cj  =  0.196854  Cg  =  0.000344 

Cg  =  0.115194  c^  =  0.019527  . 

As  before,  when  F  =  10,  we  have  P(<  lO/lO)  and  P(overcast)  =  1  -  P(<  10/ 10). 

The  error  bars  for  the  NBS  formulas  translate  to  a  maximum  computational  error 
of  1  P(:S  F/  10)  I  ^  0.  002.  In  our  example  we  set  F  =  0(1)10  in  step  2  and  show  the 
results  in  Table  4,  row  2. 

5.  OBTAINING  FREQUENCY  DISTRIBUTIONS 

In  the  example  of  Section  4,  P^  and  r  were  calculated  from  the  Revised  Uni¬ 
form  Summary  of  Surface  Weather  Observations  (RUSSWO)  sky-cover  climatology 
for  Minneapolis,  January,  1200-1400  LST.  The  procedure  for  calculating  r  is 
described  in  Section  6.7.  The  RUSSWO  provides  frequencies  for  each  tenth  of  sky 
cover  from  clear  to  overcast,  as  reported  by  the  ground  observer.  In  this  case, 
the  values  represent  the  average  of  the  three  hourly-frequencies  for  1200,  1300, 
and  1400  LST,  and  are  shown  in  Table  5,  row  1.  (See  Section  6  for  a  detailed  ex¬ 
planation  of  data  sources. )  We  assume  that  the  categories  of  tenths  as  seen  by  the 
observer  represent  equal  intervals,  such  that  the  0-tenths  category  includes  all 
individual  observations  between  0  to  0.  5  tenths,  the  1 -tenth  category  all  cases 


Table  5.  Probabilities  for  Each  Tenth  of  Cloud  Cover  for  Minneapolis, 
January  1200-1400  LST  From  (1)  RUSSWO,  (2)  the  Graphical  Method, 
and  (3)  by  Computer,  Along  With  the  rms  Errors  Between  (1),  (2)  and  (3) 


from  0.  5  to  1.5  tenths,  and  so  on,  untU,  lastly,  the  10 -tenths  (overcast)  category 
includes  all  cases  from  9.  5  to  10  tenths. 

The  graphical  and  computer  models  can  be  made  to  yield  their  results  in  the 
same  format.  First,  we  need  the  cumulative  frequencies  for  F  =  0.  5(1)9.  5-tenths 
sky-cover.  In  the  graphical  method,  they  are  obtained  by  linear  interpolation. 

For  example,  for  the  values  in  Table  4,  row  1,  P(:^  0.  5/ 10)  =  0.  5  [P(0/ 10)  + 

P(:S  1/10)]  =  0.  5  (0.  165  +  0.235)  =  0.200;  P(:£  1.  5/l0)  =  0.  5  [P(S  l/lO)  + 
P(£2/l0)]  =  0.5  (0,235  +  0.270)  =  0.253,  and  so  on.  By  computer,  the  frequen¬ 
cies  are  calculated  directly  by  setting  F  =  0.  5(1)9.  5  in  step  2  of  the  procedure 
(Section  4.2).  Since  0-tenths  coverage  for  the  observer  includes  all  cases  from 
0  to  0. 5  tenths,  we  set  the  model  result  for  P(£  0.  5/ 10)  equal  to  P(0/l0).  For 
1 -tenth  coverage,  which  for  the  observer  includes  all  cases  in  the  range  0.  5  to 
1.5  tenths,  the  model  yields  P(l/10)  =  P(sl.5/10)  -  P(£0.5/10).  We  continue 
this  process  until  we  reach  overcast,  where  P(10/l0)  =  P(>  9.  5/10)  = 

1  -  P(<9.  5/10). 

We  have  now  obtained,  via  the  model,  the  frequency  for  each  tenth  of  sky 
cover,  from  clear  to  overcast,  as  would  be  compiled  from  ground  observations. 
The  values  for  the  graphical  and  computer  methods  are  shown  in  Table  5,  rows  2 
and  3  respectively,  for  our  example.  Note  the  comparison  with  the  RUSSWO  cli¬ 
matological  data.  The  root  mean  square  (rms)  error  between  the  RUSSWO  values 
and  graphical  solution  is  0.  016,  and  between  the  RUSSWO  and  computer  solutions 
0,012. 

A  visual  comparison  for  this  example  is  shown  in  the  histogram  of  Figure 
8(b).  The  horizontal  scale  is  tenths  of  sky  cover,  and  the  vertical  scale  the  fre¬ 
quency  of  occurrence  of  each  tenth  (in  percent).  The  bars  represent  the  frequency 
distribution  of  sky  cover  as  obtained  from  the  RUSSWO,  and  the  lines  within  the 
bars  the  frequency  distribution  from  the  computer  solution.  Similar  comparisons 
are  made  for  other  RUSSWO  locations  in  Figures  8(a),  8(c),  9,  and  10.  The 
names  of  these  locations  and  times  of  observations,  as  well  as  the  RUSSWO  and 
computer  derived  cumulative  frequencies  and  rms  errors  between  them  are  given 
in  Table  6.  Figures  8  through  10  give  examples  for  different  combinations  of  r 


Table  6.  Names  of  the  Locations  and  Times  of  Observations  From  Figures  8  to  10.  Values  include 
cumulative  probabilities  of  fractional  coverage,  and  rms  errors  between  the  RUSSWO  values  and  mo 
timates.  RUSSWO  values  are  top  lines  and  model  values  bottom  lines 


rms  =  0.  014 


TENTHS  COVERAGE  TENTHS  COVERAGE  TENTHS  COVERAGE 

Figure  8.  Histograms  of  the  Frequency  Distribution  of  Total  Sky-Cover  for  Locs 
tions  With  Medium-Range  Scale -Distances.  The  bars  represent  the  RUSSWO  cli 
matologies  and  the  lines  the  computerized  model  solutions 

and  P^.  In  each  figure,  the  first  example  (a)  is  a  location  with  a  low  mean  sky- 
cover  (Pq  s  0.  3),  the  second  (b)  a  location  with  moderate  sky-cover,  and  the 
third  (c)  a  location  with  considerable  cloudiness  (P^  sO.  75).  Figure  8  contains 
cases  with  medium-range  scale  distances  (2.  0  <  r  <  2.5),  Figure  9  small 


Figure  9.  Histograms  of  the  Frequency  Distribution  of  Total  Sky-Cover  for  Loca¬ 
tions  With  Small  Scale-Distances.  The  bars  represent  the  RUSSWO  climatologies 
and  the  lines  the  computerized  model  solutions 
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Figure  10.  Histograms  of  the  Frequency  Distribution  of  Total  Sky-Cover  for  Lo¬ 
cations  With  Large  Scale -Distances,  "rfie  bars  represent  the  RUSSWO  climatolo¬ 
gies  and  the  lines  the  computerized  model  solutions 


scale -distances  (r  £0.99),  and  Figure  10  large  scale -distances  (r  s  4.00).  With¬ 
in  each  category,  the  cases  were  chosen  at  random.  The  model  fits  well  in  each 
case.  Thus,  once  the  scale  distance  and  mean  sky-cover  are  known  for  a  parti¬ 
cular  location,  the  probability  distribution  is  known  to  within  a  high  degree  of 
accuracy.  This  becomes  a  very  efficient  means  of  data  storage. 

Notice  the  effect  of  the  model  parameter,  r,  on  the  distributions.  As  r  in¬ 
creases  (going  from  Figure  9  to  Figure  10),  the  frequency  of  the  clear  and  over¬ 
cast  categories  increases  while  the  probability  of  fractional  coverage  decreases. 
The  most  striking  example  is  between  Figures  9(b)  and  10(b).  Both  locations  have 
almost  identical  mean  sky-covers  (Pq~  0.5),  and  yet  the  sky-cover  distribution 
goes  from  a  bell  shape  in  the  case  of  Figure  9(b),  where  r  =  0.  55,  to  U-shaped  in 
Figure  10(b),  where  r  =  5.  63.  This  is  consistent  with  the  definition  of  r,  which  is 
the  distance  over  which  p(s)  is  approximately  0.  99.  Thus  small  r-values  indicate 
relatively  small-scale  cloud  patterns,  such  as  summertime  cumulus,  in  contrast 
to  the  large  r-values  indicative  of  the  broad,  synoptic -scale  systems. 

The  global  maps  in  this  atlas  have  been  analyzed  on  the  synoptic  scale,  with 
mesoscale  and  local  effects  being  smoothed  out.  Therefore  it  will  not  be  possible, 
in  many  cases,  to  know  the  probability  distribution  at  any  location  on  the  atlas 
with  the  same  degree  of  accuracy  as  in  Hie  examples  of  Figures  8  through  10,  but 
it  will  be  possible  to  know  the  general  shape  of  that  distribution  within  the  region 
in  question  (e.g. ,  U-,  L-,  J-shaped)  and  see  how  that  shape  changes  from  region 
to  region.  Figure  11  depicts  the  different  distributions  for  various  mean 


22 


Figure  11.  Frequency  Distributions  in  Tenths  of  Sky-Cover,  From  Clear  to 
Overcast,  for  Various  Mean  Sky-Covers  (Pq)  and  Scale -Distances  (r/km). 

The  distributions  for  Pq  >  0.  5  are  the  mirror  images  of  the  complementary 
value  of  Pq 

sky-covers  and  scale  distances.  Distributions  for  mean  sky-covers  greater  than 
0.  5  are  the  mirror  images  of  the  complementary  means.  This  figure  can  be  used 
to  get  a  good  idea  of  the  distribution  in  a  particular  region  without  using  the 
graphical  or  computer  methods  (Section  4). 


It  is  particularly  significant  that,  in  producing  this  atlas,  it  was  not  neces¬ 
sary  to  break  the  globe  up  into  a  complicated  set  of  regions  of  rigid  cloud  clima¬ 
tologies  with  definite  boundaries,  as  many  other  atlases  do.  Therefore,  this 
atlas  is  more  realistic,  since  changes  in  the  sky-cover  distribution  take  place 
gradually  from  location  to  location. 


6.  DATA  SOURCES,  PROCESSING,  AND  ANALYSIS 

Data  were  taken  from  the  following  sources; 

1.  Revised  Uniform  Summary  of  Surface  Weather  Observations  (RUSSWO) 

2.  National  Intelligence  Survey  (NIS)  Climatic  Series 

3.  U.  S.  Navy  Marine  Climatic  (NAV)  Atlas  of  the  World 

4.  Woods  Hole  Oceanographic  Institution  (WHOI)  climatic  data 

5.  Summary  of  Synoptic  Meteorological  Observations  (SSMO) 

6.  World  Survey  of  Climatology 

6.1  RUSSWO 

The  RUSSWOs  are  weather  summaries  compiled  by  the  USAF  Environmental 
Technical  Applications  Center  and  available  through  the  National  Technical  Infor¬ 
mation  Service  (NTIS).  We  obtained  RUSSWO  sky-cover  climatologies  for  27  3 
stations.  107  in  the  U.  S.  and  166  abroad,  including  Canada.  Figure  12  is  an  ex¬ 
ample  of  the  data  format.  The  monthly  observations  are  grouped  into  eight  3 -hr 
periods,  given  in  Local  Standard  Time  (LST).  starting  with  00-02  LST  and  ending 
with  21-23  LST.  The  percent  frequency  of  total  sky-cover  is  given  in  tenths, 
from  clear  to  overcast;  each  value  represents  an  average  of  the  three  hourly-fre¬ 
quencies.  Also  shown  are  the  mean  sky-cover  (P^  in  the  BSW  model)  and  number 
of  observations  for  each  time  period,  mean  frequencies  for  the  entire  month,  and 
the  period  of  record.  The  same  format  is  used  for  all  12  months.  For  the  atlas, 
data  were  selected  for  the  mid-season  months  for  the  time  periods  00-02  .  06-08. 
12-14.  and  18-20  LST.  was  taken  directly  from  the  summaries.  The  scale 
distances  were  calculated  using  a  procedure  described  later  in  this  section.  To 
use  this  procedure,  the  data  had  to  be  converted  into  cumulative  frequencies.  To 
do  this,  we  assumed  that  the  categories  of  tenths  represent  equal  intervals  that 
cover  the  entire  spectrum  from  clear  to  overcast,  as  described  in  Section  5. 
Table  7  shows  the  cumulative  frequencies  calculated  from  Figure  12  for  January, 
00-02  LST,  for  Savannah,  Ga. 


AIR  WEATHER  SERVICE/MAC 


RUSSWO  Data  for  Savannah,  Georgia;  January 


Table  7.  Cumulative  Frequency  of  Total  Sky-Cover  for  Savannah,  Georgia, 
in  January,  for  00-02  LST  (From  Figure  12) 


X  (tenths) 

0.  5 

1.  5 

2.5  3.  5  4.  5  5.  5  6.5 

m 

8.  5  9.  5 

P(^x) 

0.349 

0.  392 

0.437  0.476  0.  509  0.526  0.552 

0.  583 

0-628  0.660 

6.2  NIS 

Section  23  of  the  NIS,  Weather  and  Climate,  was  our  largest  data  source, 
with  1846  stations  world-wide.  The  intelligence  community  updated  some  of  these 
data  until  the  early  70's,  although  for  many  stations  the  updates  ended  even  ear¬ 
lier.  The  station  data  consist  of:  a  mean  sky-cover;  the  probability  of  being  less 
than  or  equal  to  a  certain  fractional  coverage,  such  as  2/8  or  3/ 10;  and  the  prob¬ 
ability  of  exceeding  a  certain  coverage,  such  as  6/8  or  7/10.  These  values  gen¬ 
erally  appear  for  several  times  daily  for  each  of  the  twelve  months,  although 
some  countries  only  report  daily  averages.  Since  the  observing  time  varies  from 
country  to  country,  the  data  had  to  be  fit  into  the  atlas  time  periods  that  corre¬ 
spond  most  closely  to  the  time  of  the  observation.  The  only  change  that  needed  to 
be  mad«  in  these  data  for  estimating  the  scale  distances  was  to  convert  the  second 
proba  Uity  value  to  a  cumulative  frequency  by  taking  the  complement. 

6.3  NAV  Atlag 

The  NAV  atlas  (NAVAIR  50-1C-528)  provided  cloud  climatologies  for  222 
ocean  locations  around  the  world.  Most  of  the  data  are  from  ship  observations. 
Since  ships  are  always  moving,  point  statistics  were  obtained  by  collecting  a  suf¬ 
ficient  number  of  observations  in  an  area  small  enough  to  permit  an  approxima¬ 
tion.  The  data  are  given  in  the  form  of  cumulative  percent  frequency  curves 
from  0-  to  8-eights  coverage.  The  values  are  daily  averages,  and  are  given  ♦’or 
each  month.  The  cumulative  frequency  of  each  eighth  of  sky  cover  was  simply 
obtained  from  the  curves.  The  mean  sky-cover  was  then  estimated  from  the  fre¬ 
quency  distribution.  Since  ships  tend  to  avoid  bad  weather  when  possible,  these 
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data  may  be  biased  toward  good  weather  samples. 


20.  Ouayle,  R.  G.  (1974)  A  climatic  comparison  of  ocean  weather  station  and 
transient  ship  records.  Mariners  Weather  Log  18  (No.  5),  NOAA  Envi¬ 
ronmental  Data  Service. 


6.4  WHOI 


WHOI  compiled  their  mean  sky-cover  siunmaries  from  a  collection  of  surface 
marine  observations  held  by  the  National  Climatic  Center,  Asheville,  N.  C.  The 
data  were  taken  primarily  by  ships  and  within  the  established  shipping  lanes  of  the 
Atlantic  Ocean  during  the  period  1942-1972.  The  values  are  daily  and  spatial 
averages,  calculated  over  all  valid  observations  in  a  month  within  a  Marsden 
Square,  and  plotted  at  the  center  point  of  a  square.  (A  Marsden  Square  covers  an 
area  10°  latitude  by  10°  longitude.)  There  are  116  squares  partially  or  totally 
within  the  North  and  South  Atlantic  Ocean,  and  each  of  the  midseason  months  had 
mean  sky-covers  calculated  for  most  of  these  squares.  Within  any  single  10° 
square,  the  spatial  coverage  of  observations  may  vary  considerably  depending  on 
the  location  and  number  of  shipping  lanes  and  the  fraction  of  the  square  covered 
by  land  areas.  As  with  the  NAV  atlas  data,  there  may  be  a  bias  toward  good 
weather  samples. 

6.5  SSMO 

The  SSMO  is  produced  by  the  Naval  Oceanography  Command  Detachment  at 
Asheville,  N.  C.  It  is  a  series  of  coastal  marine  summaries  developed  from  ship 
observations  at  or  near  ports  world-wide.  The  mean  sky-cover  is  given  once  per 
day  for  each  month  of  the  year.  Since  the  data  are  based  upon  observations  made 
by  ships  in  passage,  the  geographic  location  of  each  value  is  the  central  position 
of  the  observations  within  an  area.  Values  were  available  from  443  locations. 

6.6  World  Survey  of  Climatology 

2 1 

The  World  Survey  of  Climatology  contains  an  appendix  with  climatic  tables 
for  125  stations  in  the  Soviet  Union.  The  sources  of  these  data  are  also  provided. 
Mean  sky-cover  is  given  three  times  daily  (0700,  1300,  and  1900  LST)  for  each 
month  at  each  station,  and  were  fit  into  the  atlas  time  periods  0600-0800,  1200- 
1400,  and  1800-2000  LST,  respectively. 

6.7  Calculation  of  Scale  Diatancea 

We  estimated  the  scale  distances  for  the  RUSSWO,  NAV  atlas,  and  NIS  data 
sets.  The  other  data  do  not  contain  probabilities  of  fractional  coverage,  which 
are  necessary  in  determining  r.  Our  procedure  was  to  use  Eq.  (8)  to  calculate 


21.  Lydolph,  Paul  E.  (1977)  Climates  of  the  Soviet  Union.  World  Survey  of  C^- 
matology  7:363-427,  Elsevier  Scientific  Publishing  Co. ,  Amsterdam,  ‘The 


2 

r  by  setting  A  =  2424  km  ,  and  z  equal  to  the  value  obtained  from  Eq.  (9).  The 

latter  equation  cannot  be  solved  to  yield  z  directly,  and  therefore  z  was  estimated 
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using  an  iterative  process  that  converges  on  z.  At  each  location  where  an  esti¬ 
mate  of  z  was  made,  the  value  was  based  on  the  geometric  average  of  several 
solutions.  For  example,  in  the  case  of  the  RUSSWO  data,  the  climatologies  con¬ 
sist  of  10  cumulative  frequency  values,  P(<  F),  where  F  =  0.  5(1)9.  5  tenths  of 
coverage,  and  a  mean  sky-cover,  P^.  The  probabilities  were  converted  into  END 
values,  P(^  F)  —  y'  and  P^  —  1  -  y^.  Ten  estimates  of  z  were  then  made  from  the 
10  pairs  of  y'  and  F  values,  and  their  geometric  mean  (z)  taken.  Outliers  were 
then  eliminated  by  omitting  z  when  |z  -  z|  >  1,  and  recalculating  z  if  necessary. 
We  then  plugged  z  into  Eq.  (8)  to  obtain  r.  The  same  procedure  was  used  to  esti¬ 
mate  r  from  the  NAV  atlas  and  NIS  data.  In  the  case  of  the  former,  there  were 
eight  pairs  of  y*  and  F  since  the  sky  cover  is  given  in  octas,  and  in  the  case  of 
the  latter,  two  pairs  of  y'  and  F. 

The  validity  of  calculating  r  using  just  two  pairs  of  y’  and  F  was  tested  by  re¬ 
calculating  r  for  selected  RUSSWO  stations  using  only  the  two  pairs.  The  new  r 
values  compared  favorably  with  those  calculated  using  the  10  pairs  of  y'  and  F; 
differences  were  within  the  noise  level.  The  conclusion  is  that  our  method  for 
calculating  r  is  as  valid  for  the  NIS  data  as  for  the  RUSSWO  and  NAV  atlas  data. 

6.8  Plotting  and  Analysis 

The  atlas  contains  a  set  of  20  global  maps  for  each  parameter,  P^  and  r,  16 
of  which  are  for  the  time  periods  0000-0200,  0600-0800,  1200-1400,  and  1800- 
2000  LST  for  the  midseason  months,  and  4  of  which  are  monthly  means,  also  for 
the  midseason  months.  Data  for  the  maps  depicting  diurnal  variations  in  P^ 
[Figures  13(a)  to  28(a))  came  from  the  RUSSWO,  World  Survey  of  Climatology, 
and  most  of  the  NIS.  The  rest  of  the  NIS,  along  with  the  NAV,  WHOI,  and  SSMO 
data  sets,  were  used  for  the  seasonal  P^  maps  [  Figures  29(a)  to  32(a)] .  Only  the 
RUSSWO,  NAV,  and  NIS  data  could  be  used  to  calculate  r,  since  they  were  the 
only  data  sets  containing  the  fractional  as  well  as  the  mean  sky-covers.  The 
RUSSWO  and  most  of  the  NIS  were  used  for  the  diurnal  depictions  of  r  [Figures 
13(b)  to  28(b)] ,  while  the  remainder  of  the  NIS  along  with  the  NAV  data  were  used 
for  the  seasonal  maps  [Figures  29(b)  to  32(b)]. 

Due  to  the  split  between  diurnal  and  seasonal  map  sets,  missing  data,  and 
regions  with  no  data,  each  map  contains  gaps  where  no  analyses  could  be  per¬ 
formed,  as  outlined  by  the  heavy  lines  in  the  analyses.  In  general,  the  diurnal 
climatologies  are  found  mostly  over  the  landrnasses,  with  the  seasonal  climatol¬ 
ogies  mostly  over  the  oceans. 


28 
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The  data  were  plotted  by  computer  on  Gringorten  polar  equal  area  maps, 
and  analyses  were  done  by  hand.  The  contours  for  are  whole  tenths  of  sky 
cover,  beginning  at  1  tenth  and  going  through  9  tenths.  The  contours  for  r  are 
as  follows:  0.25,  0.30,  0.40,  0,60,  0.80.  1.  1.  1,5,  2.0,  3.0,  6,  0,  12,0,  and 
25.  0  km.  The  reason  for  the  gradual  change  in  contour  interval  from  0.  05  to 
13  km  is  that  small  r-values  are  more  sensitive  to  changes  in  the  sky-cover  fre¬ 
quency  than  large  values.  To  illustrate  this  fact,  we  look  at  the  graphs  (Figures 
2  to  7).  The  y^  curves  have  greater  slopes,  in  general,  for  large  z-values  (small 
r-values)  than  for  small  z-values  (large  r-values).  In  addition,  an  interval  on  the 
z  scale  at  the  high  end  corresponds  to  a  smaller  change  in  the  values  of  r  than  an 
Identical  interval  on  the  low  end  of  the  z  scale  [  see  Eq.  (8)) .  For  instance,  from 
z  =  8  to  z  =  7,  r  goes  from  0.2  to  0.4  km,  whereas  from  z  =  3toz  =  2,  r  goes 
from  6  to  12  km.  In  both  instances  ^z  =  1,  but  ^r  =  0.  2  km  in  the  first  case  and 
6  km  in  the  second.  Thus,  with  a  horizontal  axis  uniform  in  r  instead  of  z,  the 
change  in  the  slopes  of  the  y^  curves  would  be  greatly  accentuated.  In  order  to 
compensate  for  this  effect,  the  contour  intervals  for  r  were  designed  to  permit 
more  precise  map  readings  for  small  values  than  for  large  values. 


22.  Gringorten,  I.  I.  (1981)  Mapping  the  Climate,  AFGL-TR-81-0015, 
AD  A102904. 
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Figure  23b.  Analysis  of  Scale  Distance  (r)  for  July,  1200  to  1400  LST 
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